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Abstract. The isotropy of gamma-ray bursts collected 
in Current BATSE catalog is studied. It is shown that 
the quadrupole term being proportional to ~ sin 2b sin I is 
non-zero with a probability 99.9%. The occurrence of this 
anisotropy term is then confirmed by the binomial test 
even with the probability 99.97 %. Hence, the sky dis- 
tribution of all known gamma-ray bursts is anisotropic. 
It is also argued that this anisotropy cannot be caused 
exclusively by instrumental effects due to the nonuniform 
sky exposure of BATSE instrument. Separating the GRBs 
into short and long subclasses, it is shown that the short 
ones are distributed anisotropically, but the long ones 
seem to be distributed still isotropically. The character of 
anisotropy suggests that the cosmological origin of short 
GRBs further holds, and there is no evidence for their 
Galactical origin. 
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1. Introduction 

The true physical nature of gamma-ray bursts (GRBs) is 
one of the tantalizing enigmas of the recent astrophysics. 
Although since their first detection (Klebesadel et al. 
1973) there were several suggestions trying to give a clear 
explanation for their origin, no definite answer has been 
given yet (cf. Paczyhski 1995). Recently, the successful 
identifications made by the Beppo-SAX satellite, followed 
by the detection of optical counterparts (van Paradijs et 
al. 1997), seem to give a firm support for the models aim- 
ing to explain the bursts by merging neutron stars (Usov 
& Chibisov 1975; Rees & Meszaros P. 1994; Meszaros P. 
& Rees 1997) and seem to put them definitely into the 
cosmological distances. The alternative Galactical origin 
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seems to be ruled out (a survey of the question of dis- 
tances may be found, e.g., in Paczyhski 1995). However, 
the small number of optically identified events is far from 
being enough to characterize the properties of the whole 
burst population. On the other hand, the existence of cos- 
mological distances of GRBs seems to be definite. 

In addition, even before this identification, indirect ob- 
servational evidences were known for the cosmological ori- 
gin. These evidences were based mainly on the modifica- 
tions of < V/Vmax > test (cf. Norris et al. 1994; Meszaros 
P. & Meszaros A. 1995; Norris et al. 1995; Nemiroff 1995; 
Horvath et al. 1996), and on the study of the time di- 
latation (cf. Norris et al. 1995; Meszaros A. & Meszaros 
P. 1996; Stern 1996; Meszaros A. et al. 1996; Che et al. 
1997). A further important indirect support of cosmolog- 
ical origin is based on the observed isotropy on the sky 
(Briggs 1993; Syer & Saha 1994; Briggs 1995; Tegmark et 
al. 1996a, b; Briggs et al. 1996). All these papers suggest 
that the angular distribution is isotropic, because there are 
no statistically significant departures from the isotropy. 
Also the separations of GRBs into the different subclasses 
either due to the duration (Kouveliotou et al. 1993; Belli 
1995; Dezalay et al. 1996) or due to the fluence on channel 
above the 300 keV (Pendleton et al. 1997) do not change 
the situation; the proposed subclasses alone also seem to 
be distributed isotropically. Hence, it can well occur that 
the total number of observed GRBs is a mixture of a wide 
variety of physically different objects, but all GRBs should 
be at cosmological distances due to their isotropic angular 
distribution on sky, and due to other direct and indirect 
supports. 

In this paper we find a clear anisotropy of all GRBs, 
and then separately of the short GRBs, too. On the other 
hand, the long GRBs seem to be still distributed isotrop- 
ically. To do this we will use both the standard analysis 
based on the spherical harmonics and also the so called 
binomial test. It may seem that all this is an argument 
against the cosmological origin of short GRBs. It is shown 
that this is not the case. 
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The paper is organized as follows. In Section 2 the 
key ideas of the used statistical tests are recapitulated. 
In Section 3 the anisotropy of all GRBs is demonstrated. 
In Section 4 the anisotropy of short and the isotropy of 
long GRBs is shown. Section 5 discusses the results, and 
argues for the cosmological origin of short GRBs. Finally, 
in Section 6 there is a summation of results. 

2. Mathematical skeleton of the problem 

Testing the isotropy on the celestial sphere one may use 
several methods. Nevertheless, strictly from the mathe- 
matical point of view, the necessary condition for the 
isotropy is the stochastic independency of the sky distri- 
bution of the bursts on their observed physical properties. 
It means that, if f(b, I, X\, x n )dFdxi...dx2 is the proba- 
bility of finding an object in the dF = cos 6 dl db infinites- 
imal solid angle and in the (x\,x\ + dx\, x n ,x n + dx n ) 
interval, one must have 

f(l,b,x u ...,x n ) = u{l,b)g{x u ...,x n ). (1) 

Here < I < 360°, -90° < b < 90° give the celestial 
positions in Galactical coordinates, x n (n > 1) measure 
the physical properties (peak fluxes, fluences, durations, 
etc..) of GRBs and g is their probability density. It may 
well be assumed that in the case of isotropy the distribu- 
tion of GRBs fulfils this equation (cf. Briggs et al. 1996; 
Tegmark et al. 1996a, b). 

However, statement (1) is only a necessary but not a 
sufficient condition for isotropy. Isotropy means that also 
oj(1 7 b) = l/(47r). Hence, in this case, for N observed GRBs 
the events dN — NuidF, i.e. the expected number of GRBs 
in an infinitesimal solid angle, is not depending on [l,b]. 
In other words, the isotropy means that the probability of 
observation of a burst in a solid angle < f2 < 4w (fi is in 
steradians) is given by n/(4ir) and is independent on its 
location on the celestial sphere. This follows immediately 
from (1), if one does integration over I and b to obtain, 
first, the solid angle O, and, second, the whole sky. Then 
the ratio of two results gives fi/(47r), and the concrete 
form of g is unimportant. 

The most frequently used procedure to test the 
isotropy of GRBs is based on the spherical harmonics 
(Briggs 1993, 1995; Briggs et al. 1996; Tegmark et al. 
1996a, b). The key idea is the following. In general case one 
may decompose the function iv(b, I) into the well-known 
spherical harmonics. One has: 

oj(b,l) = (4tt)- 1 / 2 wo- 
(3/(47r)) 1 / 2 (wi,_i cos6sinZ — u>i,i cos6cosZ + w^o sin6)+ 
(15/(167r)) 1,/2 (aj2,-2 cos 2 6sin2Z + w 2 ,2 cos 2 6cos2Z— 
u)2,-i sin 26 sin I — w 2 ,i sin26cos/)+ 

(5/(167r)) 1 / 2 tj2,o(3 sin 2 b—l)+higher order harmonics. (2) 



The first term on the right-hand side is the monopole 
term, the following three ones are the dipole terms, the fol- 
lowing five ones are the quadrupole terms (cf. Press et al. 
1992; Chapt. 6.8). Nevertheless, u is constant for isotropic 
distribution, and hence on the right-hand side any terms, 
except for ojq, should be identically zeros. To test this hy- 
pothesis one may proceed, e.g., as follows. Let there are 
observed N GRBs with their measured positions [bj,lj] 
(j = 1, 2, N). In this case w is given as a set of points 
on the celestial sphere. Because the spherical harmonics 
are orthogonal functions, to calculate the u>q coefficients 
one has to compute the functional scalar products. For 
example, w 2 ,-i is given by 

/ 15 \V2 r r 

U2,-i = -[j£-) / cosbdb / ui(l, b) sin 26 sin Z dl 

-tt/2 
N 

= -(15/(167r)) 1 / 2 Ar- 1 ^sin26 i sinZ J . (3) 

Because u> is given only in discrete points, the integral 
is transformed into an ordinary summation (cf. Kendall 
& Stuart 1969, p. 16). In the case of isotropy one has 

N 

LU2.-1 = 0, and hence N^ 1 sin26j sinZj = 0. Therefore, 

i=i 

the expected mean of SU126J sin^- values is zero. One has 
to proceed similarly to any other coefficient. 

In order to test the zero value of, e.g., cj 2 .-i one has 
to calculate, first, sin26 J - sinij for any j = 1, 2, N, and, 
second, the mean, standard deviation and Student's t vari- 
able (cf. Press et al. 1992, Chapt. 14). Finally, third, 
one has to ensure the validity of zero mean from Stu- 
dent test. As far as it is known, no statistically significant 
anisotropics of GRBs were detected yet by this procedure 
(cf. Briggs et al. 1996; Tegmark et al. 1996a, b). 

Nevertheless, there are also other ways to test the 
isotropy. An extremely simple method uses the binomial 
distribution. In the remaining part of this section we ex- 
plain this test (see also Meszaros A. 1997). 

In order to test the anisotropy by this method one may 
proceed as follows. Let us given an area on the sky defined 
by a solid angle < ft < An (Q is in steradians). In the 
case of isotropy the p probability to observe a burst within 
this area is p = fl/ (47r). Then, obviously, q = 1 — p is the 
probability to have it outside. Observing > bursts 
on the whole sky the probability to have k bursts (it may 
be k — 0, 1, 2, N) within O is given by the binomial 
(Bernoulli) distribution taking the form 

This distribution is one of the standard probability dis- 
tributions discussed widely in statistical text-books (e.g. 
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Trumplcr & Weaver 1953; Kendall & Stuart 1969, p.120; 
about its use in astronomy see, e.g., Meszaros 1997). The 
expected mean is Np and the expected variance is Npq. 
One may also calculate the integral (full) probability, too, 
by a simple summation. 

In our case we will consider N GRBs, and we will test 
the hypothesis whether they are distributed isotropically 
on the sky. Assume that k b s is the observed number of 
GRBs at the solid angle £1. If the apriori assumption is 
the isotropy, i.e p = Q/(47r), then one may test whether 
the observed number k b s is compatible with this apriori 
assumption. Of course, any < k b s < N can occur with 
a certain probability given by the binomial distribution. 
But, if this probability is too small, one hesitates seriously 
to accept the apriori assumption. 

Consider the value | k b. s . — Np \= k a . The value k a 
characterizes "the departure" of k bs from the mean Np. 
Then one may introduce the probability 

P(N,k obs ) = 

1 - P P MN, (Np + ho)) + P P ,tot(N, (Np - ho)). (5) 

P(N,k b s ) is the probability that the departure of k a b s 
from the Np mean is still given by a chance. 

In order to test the isotropy of the GRBs celestial dis- 
tribution we will divide the sky into two equal areas, i.e 
we will choose p = 0.5. It is essential to note here that nei- 
ther of these two parts must be simply connected compact 
regions. 

3. Anisotropy of all GRBs 

In order to test the isotropy of GRBs we will test the three 
dipole and five quadrupole terms in accordance with the 
method described in the previous Section. We will con- 
sider all GRBs that have well-defined angular positions. 
Up to the end of year 1997 there were 2025 such objects 
at Current BATSE catalog (Meegan ct al. 1997; Paciesas 
et al. 1998). The results are summarized in Table 1. 

We see that, except for the terms defined by u>2,-i and 
W2,-2, the remaining six terms may still be taken to be 
zero. This means that there is a clear anisotropy defined by 
term ~ sin 2b sin I. The probability that this term is zero 
is smaller than 0.27%. It is practically sure that the sec- 
ond quadrupole term being proportional to ~ cos 2 b sin 21 
is non-zero, too. Nevertheless, we will not deal with this 
second term in this paper, because the purpose of this pa- 
per is to demonstrate only qualitatively the existence of 
anisotropies. 

The anisotropy defined by o> 2 ,-2 ma y be defined in an- 
other form, too. This quadrupole term has a positive sign, 
when both b and I have the same signs, and has a neg- 
ative sign, when I and b have opposite signs. Therefore, 
let us define two parts of the sky having the same sizes 





t 


% 


Wl,-1 


1.51 


13.4 


Wl,l 


1.77 


7.7 


Wl,0 


0.71 


47.7 


W2,-2 


2.76 


0.6 


<^2,2 


1.54 


12.1 


W2.-1 


3.26 


0.1 


W2,l 


0.98 


33.3 


W2,0 


0.36 


71.9 



Table 1. Student test of the dipole and quadrupole terms of 
2025 GRBs. In the first column the coefficients denned in Eq.2 
are given. In the second column the Student t is provided. The 
third column shows the probability that the considered terms 
are still zeros. 



(27T steradians). The first one is defined by Galactical co- 
ordinates b > 0, I > 180° and b < 0, < I < 180°. This 
means that this first part is in fact composed from two 
separated " sky-quarters" . The second part is then given 
by b > 0, < I < 180° and b < 0, I > 180°. This means 
that this second part is also given by two separated "sky- 
quarters" . Then the detected quadrupole anisotropy sug- 
gests that there should be an essential difference, e.g., in 
the number of GRBs in these two parts. 

In order to test again this expectation we will do the 
test based on Bernoulli distribution. We divide the whole 
sky into these two parts, and hence we expect a Bernoulli 
distribution with p = 0.5 for N = 2025. As it is noted 
in Section 2, it is certainly allowable that these parts are 
composed from several subregions. 

A straightforward counting of GRBs in these regions 
shows that 930 GRBs are in the first one and 1095 are in 
the second. (Note that no GRBs had coordinates exactly 
either b = 0, or b = ±90°, or / = 0, or I = 180°. This means 
that in the paper no problems have arisen from the fact 
that the boundaries of "sky-quarters" were not taken into 
acount.) Assuming p — 0.5 the binomial (Bernoulli) test 
gives a 0.03% probability that this distribution is caused 
only by a chance. Hence, the relatively smaller number 
in the first region compared with the second one is not 
a chance, and the distribution of all GRBs is anisotropic 
with a certainty. 

Clearly, concerning the consequences of the intrinsic 
anisotropy of GRBs, one must be still careful. Instru- 
mental effects of the BATSE experiment may also play 
a role, and in principle it can also occur that the de- 
tected anisotropy is caused exclusively by instrumental 
effects. To be as correct as possible, one may claim that, 
in essence, there can be three different causes of this ob- 
served anisotropy: a. The anisotropy is purely caused by 
the BATSE's nonuniform sky exposure (in other words, 
the intrinsic angular distribution of GRBs is still isotropic, 
and the observed anisotropy is a pure instrumental ef- 
fect); b. The anisotropy is purely given by the intrin- 
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sic anisotropy of GRBs, and the instrumental effects are 
unimportant; c. The anisotropy is given both by instru- 
mental effects and also by the intrinsic anisotropy. To be 
sure that there is also an intrinsic anisotropy of GRBs, 
one must be sure that the possibility a. does not occur. In 
what follows, when we will speak about "the possibility 
a." , we will consider this one. 

It is well-known that the sky exposure of the BATSE 
instrument is nonuniform. This question is described and 
discussed in several papers (cf. Briggs 1993; Fishman et 
al. 1994; Tegmark et al. 1996a,b; Briggs et al. 1996). The 
BATSE sky coverage depends on the declination only in 
the equatorial coordinate system (Tegmark et al. 1996b) 
in such a manner that the probability of detection is about 
10% higher near the pole than near the equator. This be- 
haviour in Galactic coordinates predicts excess numbers 
of GRBs just in sky-quarters given by b > 0, < I < 180° 
and b < 0, I > 180°. Hence, in principle, it is well possible 
that the observed anisotropy is caused by a pure instru- 
mental effect. The purpose of Section 4 is to show that this 
is not the case, and the possibility a. should be excluded. 

4. Different distribution of short and long GRBs 

If there is an intrinsic isotropy indeed, then Eq.(l) will be 
further fulfilled, and ui(b,l) itself will reflects the nonuni- 
formity of sky-coverage. Then at the first part the number 
of observed GRBs should be smaller, because the integra- 
tions of uj(1, b)dF giving the first and second part, respec- 
tively, do not give the same values. Their ratio should be 
~ 930/1095 = 0.85. This ratio should be obtained for any 
subset of GRBs, when the choice of this subset is based on 
some physical properties of the bursts, because the func- 
tion g does not enter into the calculation for X. In other 
words, if the level of anisotropy depends on the physical 
parameters of GRBs, then "the possibility a." should be 
excluded. 

To be extra cautious it is also necessary to remember 
that some bias may arise for the dimmest GRBs, because 
for them it is not necessarily true that the total exposure 
time is exactly proportional to the observed number of 
sources due to the varying threshold limit of BATSE (see 
Tegmark et al. 1996a for the discussion of this question). 
To avoid this bias the simplest procedure is not to take 
into account the dimmest GRBs. 

All this allows to exclude the possibility a. quite sim- 
ply. To do this it is necessary to take some subsets of 
GRBs, and to verify for any of them that the ratio X is 
roughly 0.85. For security, it is also necessary to omit the 
dimmest GRBs. 

For the sake of maximal correctness we will not use 
any ad hoc criterions to define such subsets, but we will 
exclusively use criterions which were introduced earlier by 
others. First, we exclude any GRBs having the peak fluxes 
on 256 ms trigger smaller than 0.65 photons/ (cm 2 s). The 
truncation with this threshold is proposed and used in 



Pendleton et al. (1997). Second, from the remaining GRBs 
we exclude GRBs which have no defined duration Tgo (for 
the definition of this duration see Kouveliotou et al. 1993). 
Third, we also exclude the bursts which have no /3 value, 
which is the fluence on the energy channel [100, 300] keV. 
These truncations are also necessary, because we will con- 
sider the subsets, for which the criterions use T go and /3. 

932 GRBs remain, and from them 430 are at the first 
part, and 502 are at the second one. Hence, here X = 
430/502 = 0.86. There is no doubt that this "truncated" 
sample of GRBs is distributed similarly to that of whole 
sample with 2025 GRBs. The probability that this dis- 
tribution is given by a chance is here 2%. (Note that, of 
course, here the same X does not give the same probabil- 
ity, because there is a smaller number in the sample.) 

932 GRBs will be separated, first, into the "short" and 
"long" subclasses (cf. Kouveliotou et al. 1993), and, sec- 
ond, into the No-High-Energy (NHE) and High-Energy 
(HE) bursts (Pendleton et al. 1997). 

The boundary between the short and long bursts is 
usually taken for T 90 = 2s (Kouveliotou et al. 1993; Belli 
1995; Dezalay et al. 1997). Nevertheless, this boundary at 
Tgo = 2 s is not so precise (e.g. in Katz & Canel (1996) 
T 90 = 10 s is used). In addition, this boundary gives no 
definite strict separation, because at class T 90 < 2 s long 
bursts, and at class Tgo > 2 s short bursts are also pos- 
sible, respectively (cf. Belli 1995). Therefore, in order to 
test more safely the distribution of two subclasses, we will 
consider also the case when the boundary is at Tgo = 10 
s. We will consider also the subsamples of GRBs having 
T 90 < 1 s, and T 90 > 15 s, respectively, because then they 
contain surely only short and long bursts, respectively. 
The results are shown in Table 2. 



sample 


N 


kobs 


(N -k oba ) 


% 


all GRBs 


932 


430 


502 


2.0 


T 90 < Is 


206 


82 


124 


0.43 


T 90 < 2s 


251 


103 


148 


0.55 


T 90 < 10s 


372 


154 


218 


0.11 


T 90 > 2s 


681 


327 


354 


32 


T 90 > 10s 


560 


276 


284 


77 


Too > 15s 


507 


247 


260 


59 



Table 2. Results of the binomial test of subsamples of GRBs 
with different durations. N is the number of GRBs at the given 
subsample, k b s is the observed number GRBs at the first part 
in this subsample, and % is the probability in percentages that 
the assumption of isotropy is still valid. 



Table 2 shows, e.g., that there are 251 GRBs with 
T 90 < 2 s, and 681 GRBs with T 90 > 2 s. Then, from 
the short GRBs 103 are at the first part of sky and 148 at 
the second one. This gives X = 103/148 = 0.70. It seems 
that at the first part there is even a smaller portion of 
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shorter GRBs than that of the all GRBs. The probability 
that this is a chance is given by 0.55%. On the other hand, 
from the long GRBs 327 are at the first part and 354 at 
the second one. This gives X = 327/354 = 0.92. Hence, 
it seems immediately that for the long GRBs the isotropy 
is still an acceptable assumption. The binomial test quan- 
tifies: there is a 32% probability that this distribution is 
given by a chance. Doubtlessly, the long GRBs are dis- 
tributed more isotropically than the short ones; there is 
no statistically significant departure from isotropy for long 
GRBs. The subsamples Tgo < 1 s and Tgo > 15 s con- 
firm this expectation; the boundary at Tgo = 10 s also 
does not change the conclusion. One may claim that the 
anisotropy of short GRBs is statistically significant, but 
for long GRBs it is not. 

Doubtlessly, the short and long subclasses are dis- 
tributed differently This also excludes the possibility a.; 
the observed anisotropy of all GRBs cannot be caused 
exclusively by instrumental effects. It is difficult to imag- 
ine an instrumental effect which leads to isotropy of long 
GRBs and to anisotropy of short GRBs. 

Pendleton et al. (1997) introduces the subclasses of 
HE and NHE bursts. The criterion depends on the ratio 
/4//3, where /3 is the fluence on energy channel [100, 300] 
keV and /4 is the fluence on the energy channel > 300 
keV. (From this it is also clear, why we needed non-zero 
/3. On the other hand, /4 can be vanishing; these GRBs 
are simply NHE bursts.) Application of this criterion is not 
so simple, because for a great portion of GRBs there are 
large uncertainties of the values of /4 due to their errors. 
Concretely, for 693 GRBs the value of /4 is bigger than 
the corresponding error of this /4; for 131 GRBs there is 
no /4; for the remaining 108 GRBs there are some values 
of /4 at the Current BATSE catalog (Meegan et al. 1997), 
but they are smaller than their errors. Hence, 693 GRBs 
can be taken as HE bursts (HE1 subsample). 131 GRBs 
having no /4 may be taken as NHE bursts (NHE1 sub- 
sample). We did binomial tests for these two subclasses. 
Separation of the remaining 108 GRBs into the HE and 
NHE subclasses is not so clear. We consider artificially the 
boundary as follows: If the value of /4 is bigger than the 
half of error, then we have HE; otherwise NHE. Applying 
this criterion we will have 168 NHE (NHE2 subsample) 
and 764 HE bursts (HE2 subsample). For them the bi- 
nomial tests were also done. The results are collected in 
Table 3. 

Table 3. gives an ambiguous result. Due to the smaller 
number in subclasses no anisotropy is confirmed yet on a 
satisfactorily high level of significance. In addition, con- 
trary to the short-long separation, there is no obvious dif- 
ference between HE and NHE classes. Hence, there is no 
obvious and unambiguous result here. 



e ri 777 r> / p 

tDlAil 1 LjJl C- 


N 


k t. 

™obs 


CAT — k t, ) 


% 


all GRBs 


932 


430 


502 


2.0 


NHE1 


131 


52 


69 


14.9 


HE1 


693 


327 


366 


14.6 


NHE2 


168 


69 


99 


2.5 


HE2 


764 


361 


403 


13.8 



Table 3. Results of the binomial test of NHE-HE subsam- 
ples of GRBs. The subsamples NHE1, HE1, NHE2, HE2 are 
explained in the text. 



5. Discussion 

The quadrupole anisotropy reported in this paper is an 
unexpected and new result. As far as it is known no 
anisotropy terms were detected yet (cf. Tegmark et al. 
1996a, b). Probably this situation was given by the fact 
that the majority of these isotropy studies concentrated 
the effort into dipole and LO2.0 quadrupole terms, which 
are expected to differ from zero, if the GRBs are arisen in 
the Galaxy. 

The essentially different angular distribution of short 
and long GRBs suggest that their separation into these 
two subclasses has a deeper cause. It is well-known that in 
average the short bursts have higher hardnesses (hardness 
= /3//2, where /2 is the fluence on energy channel [50, 
100] keV). Katz & Canel (1996) have also shown that the 
< V/Vmax > values are different; the smaller value for 
longer GRBs suggests that they are on average at bigger 
cosmological distances. Keeping all this in mind it seems 
to be definite that these two types are physically different 
objects at different cosmological scales. 

Contrary to this, we did not find any significant dif- 
ference in the angular distribution of HE and NHE sub- 
classes. This is an unclear result, because in Pendleton et 
al. (1997) it is clearly stated that the < V/V max > values 
are different for the HE and NHE subclasses, and hence 
they should also be at different distances. The isotropy 
tests do not confirm this expectation. This also means 
that the separation based on the most energetic channel 
remains unclear. In fact, the question of fourth channel is 
highly topical recently, because Bagoly et al. (1998) shows 
- independently on Pendleton et al. (1997) - that /4 alone 
is an important quantity. The question of most energetic 
channel trivially needs further study, and is planned to be 
done. 

At the end of Section 3 we pointed out that the de- 
pendence of the BATSE detection probability on the dec- 
lination may mimic some sort of anisotropy. Therefore, 
without a detailed correction in accordance with the sky 
exposure function one may state only the presence of an 
intrinsic anisotropy from the different behaviour of short 
and long bursts. The requirement of the study of this cor- 
rection for both subclasses is a triviality and is planned 
to be done in the near future. 
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One may also speculate that the short GRBs can arise 
in the Galaxy and the long ones at cosmological distances. 
(About the cosmological origin of long GRBs there seems 
to exist no doubt; see, cf., Meszaros A. et al. 1996). The 
existence of non-zero LU2.-1 and probably also of LO1.-2 
quadrupole terms with the simultaneous zeros for other 
dipole and quadrupole terms is a strange behavior for any 
objects arising in the Galaxy. E.g., it is highly complicated 
to have an Ct>2,-i term, and simultaneously not to have the 
^2,o term, if the sources have arisen in the Galaxy. Sim- 
ply, any objects in the Galaxy should have fully different 
anisotropy terms (for further details and for the survey 
of earlier studies of isotropy see, e.g., Briggs 1993, 1995; 
Briggs et al. 1996; Tegmark et al. 1996a, b; Meegan et al. 
1996). 

Remark also the following. The so called transition 
scale to homogeneity (cf. Meszaros A. 1997) is minimally 
of size ~ 300ft. -1 Mpc (h is the Hubble constant in units 
100 km/(s Mpc)). This means that up to this distance 
an inhomogeneous and anisotropic spatial distribution is 
not only possible but it is even expected. In addition, 
at the last time several observational implications both 
from the distribution of galaxies and from the anisotropy 
of cosmic microwave background radiation highly query 
the fulfilment of homogeneity and isotropy even up to 
the the Hubble scale (Lauer & Postman 1994; Slechta & 
Meszaros A. 1997; Meszaros A. k Vanysek 1997; Coles 
1998; Sylos-Labini et al. 1998). Hence, as anisotropy con- 
cernes we think that the different distribution of short 
and long GRBs and their cosmological origin are not in 
contradiction, and we further mean that all GRBs are at 
cosmological distances. 

6. Conclusion 

The results of paper may be collected as follows: 

A. The distribution of 2025 GRBs is anisotropic. 

B. This anisotropy is not caused exclusively by instru- 
mental effects. 

C. Separating GRBs into the short and long subclasses 
it is shown that the short ones are distributed anisotrop- 
ically, but the long ones can still be distributed isotropi- 
cally. 

D. No such difference was found for the HE-NHE sep- 
aration. 

E. We conjectured that the anisotropic distribution of 
short GRBs does not query their cosmological origin. 
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